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During this period, three papers were presented 1 ’ 3 (copies enclosed) and three abstracts 
were accepted for the 97 AIAA Sciences Meeting and Exhibit to be held at Reno in January, 1997. 

Progress is being made on a number of fronts. We have developed a k-£ (Enstrophy) 
model capable of predicting turbulence separation and applied it to two airfoils: NACA 0012 and 
RAE 2822 at various angles of attack and Mach numbers. Moreover, a two-equation k-£ model 
with a tensor eddy viscosity has been developed. The goal of this work is to eliminate the 
empirical damping factor, fp., that is used in defining the eddy viscosity in wall bounded flows. 
Plans call for this model to be used in calculating three dimensional turbulent flows. 

The work on compressible turbulence is progressing well. Current efforts are devoted to 
analyzing crossing shocks interaction with turbulent boundary layers. 

Finally, we have developed an alternative to the e 11 method using models similar to those 
employed in turbulence modeling. The method was implemented into boundary layer codes and is 
being incorporated into a Navier-Stokes solver. As a result of this method, the onset of transition 
and the rest of the flowfield can be obtained from a single Navier-Stokes code. 
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Abstract 

In an attempt to model a recent set of experiments 
by Skare and Krogstad dealing with equilibrium bound- 
ary layers near separation, it became clear that the tradi- 
tional &-to model was unable to reproduce the data. 
Careful investigation of the data suggested that the dif- 
fusion term in the k - equation behaves differently in the 
presence and absence of unfavorable pressure gradients. 
This led to the conclusion that an important diffusion 
mechanism was missing from the k - equation and a new 
diffusion model was proposed. The results of the new 
theory indicate marked improvement when compared 
with the data for strong adverse pressure gradient flows. 
Moreover, the new theory reproduces Stratford’s limit 
of velocity in a region of vanishing shear stress. How- 
ever, because subsonic and supersonic boundary layers 
act differently in the presence of an adverse pressure 
gradient, no adjustments in the diffusion term is recom- 
mended for supersonic flows in the presence of mild 
adverse pressure gradients. 

Introduction 

It is a known fact that traditional turbulence models 
are somewhat inadequate when it comes to predicting 
flows with an adverse pressure gradient 1 . This is a result 
of prevailing modeling practices. Whether we deal with 
one-equation, two-equation, or a stress model, current 
modeling does not address situations where adverse 
pressure gradients play a role. Thus, we use results from 
the decay of homogeneous turbulence, the log-law 
region, and asymptotic expansions along turbulent-non- 
turbuient boundaries to determine model constants 1 . 
There is no reason to expect that such practices will pro- 
duce turbulence models that are valid in the presence of 
unfavorable pressure gradients. 
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The situation has changed recently with the avail- 
ability of well documented incompressible experiments 
describing equilibrium boundary layers near 
separation 2 ’ 3 . As is shown in Ref. 3, the diffusion term 
in the k - equation behaves differently in the presence 
and absence of pressure gradients. This led us to the 
conclusion that an important diffusion mechanism is 
missing from the ik - equation. 

The objective of this investigation is to show that 
there is indeed a missing term resulting from the tradi- 
tional neglect of the term p r u? , where p' and u[ are 
the pressure and velocity fluctuations. Moreover, it 
results in improved agreement between theory and 
experiment for low speed flows. 

However, subsonic and supersonic boundary layers 
are known 4 to behave differently in the presence of 
adverse pressure gradients. As a result, a modification 
like the one suggested here is not expected to apply for 
all Mach numbers. It has been suggested that compres- 
sion associated with an adverse pressure gradient is 
associated with an increase in skin friction and a 
decrease in boundary layer thickness at supersonic 
speeds. This would suggest that some modification of 
the proposed turbulence model is required. While exam- 
ining the experiment of Fernando and Smits 5 , we dis- 
covered that a traditional k- co model was quite 
satisfactory. As a result, it is concluded that adjustment 
in the diffusion term is only required for subsonic flows. 

Analysis 

The application of an adverse pressure gradient on a 
wall bounded flow promotes separation of the boundary 
layer. Thus, we use the theoretical behavior of turbulent 
boundary layers near separation to investigate current 
modeling practice. 

The only term that requires modeling in the k - 
equation is the diffusion term 1 . Traditionally, it is mod- 
eled as 

l 

+/>«; = ( 1 ) 
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where p is the density, k is the turbulent kinetic energy, 
\h is the turbulent viscosity, and a* is a model constant. 

The applicability of this model near separation is 
evaluated for an incompressible flow. In the near wall 
region, the Wilcox k - co model reduces to 


0 = 


dP dx 
dx dy 


( 2 ) 


i l 

v , = k/'V (in 

When Eqs. (9). (10), and (11) are substituted into 
Eqs. (3) and (4), the system yields two equations in two 
unknowns. Solution of these equations reveals the 
asymptotic behavior in the region near separation. The 
general procedure entails setting 



where the turbulent viscosity is defined as 



P is the pressure, x is the shear stress, k is the turbulent 
kinetic energy, and to is the specific dissipation rate. The 
values of the model constants are 


k = Ay” , CO = By n (12) 

where A and B are functions of P' . Eq. (3) becomes 

( i \ i m-i i - - 

0 = m\m + ~ycs*K 0 F y 2 + JLp 2 y 2 ( 13 ) 

If Eq. (13) is independent of y, then 
I 1 

m ~2 = 2 = m + n ( 14 ) 
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Integration of Eq. (2) gives 


( 7 ) 
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^ dP 
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where x w is the wall shear stress. In regions where 
x w 0 » Stratford 6 showed that U has the form 


u = 2{JP r y)/K o , P‘ = (9) 

where K 0 is typically a constant of about 0.5. Equation 
(9) was later confirmed by Townsend^ using the data of 
Schubauer and Klebanoff®. 

The velocity gradient given by Eq. (9) is 


resulting in 


m = I, n = ~ 

2 

Thus. Eqs. (5). (11), (12), and (15) give 



1 

2 


Substituting above into Eq. (4) yields 


t aP Pfi 

fQy 


which simplifies to 



( 15 ) 


( 16 ) 


( 17 ) 


( 18 ) 


du _ j_r py 

dy ~ K„ l y ) 


( 10 ) 


Applying Eqs. (9) and (10) to Eq. (8) in the limit of van- 
ishing shear stress yields 


which gives 


A 



( 19 ) 


When Eqs. (15), (18), and (19) are applied to Eq. (13), 
one finds 
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same can be shown for the k - e model. 



JL a P 
P 


Noting that 


( 20 ) 


a 


I! 

P 


( 21 ) 


from the values of the constants, Eq. (20) reduces to 



( 22 ) 


Results and nisr.i«in n 

A boundary layer code 10 is employed in carrying 
out the computations. A standard k - co model and a 
modified k - to model were implemented. A k - to model 
and not a k - e was selected because it has been found 11 
that af-w model performs better in the presence of 
adverse pressure gradients. Detailed data on the condi- 
tions of Ref. 2 and Ref. 5 were used to test the turbu- 
lence models. These experiments will be denoted 
Experiment (1) and Experiment (2), respectively. 


In view of Eq. (13), the above result is a contradiction. 
This shows that the current modeling of diffusion is 
inadequate in a region near separation. 

To remedy the above situation, we supplement dif- 
fusion with an additional term. The proposed modeling 
of the diffusion term is 


1 . , 

2 f >u i u i u j + P'“j 


= f 

L dx i p p dx j\ 


where 


(23) 


h - T '/ . 2 X — , — , 

iJ ~ pi + T V = -P“7^ ( 24 ) 

and o p is a model constant. This model constant is 
determined from the same analysis performed above. 
Using the standard k - to notation, it can be shown that 



for K 0 =0.51. The value of K 0 was found by Townsend 7 
based on the data of Schubauer and Klebanoff 8 . Note 
that o p is negative, so the new term is a negative diffu- 
sion term in the near wall region for an adverse pressure 
gradient. It should be pointed out that the Stratford limit 
is valid in the near-wail region where x w —>0. The flow 
under considerauon is attached. Thus the value pre- 
sented in Eq. (25) is to be treated as a crude estimate of 
<J p . After computadons were made for a strong adverse 
pressure gradient flow, we set o' = -5.13 to best match 
the experiment. 

It should be noted that Saffman' attempted to repro- 
duce the above limit using his k - co 2 model. He found 
that in order for his model to reproduce Eq. (9), the vari- 
ous flow parameters must take on imaginary values. The 


Experiment m 

Experiment (1) reports an equilibrium boundary 
layer in a strong adverse pressure gradient flow at a 
Reynolds number per length of approximately 1.2 mil- 
lion. An equilibrium flow exhibits similarity in the outer 
region of the boundary layer 7 . The equilibrium region is 
maintained between x c = 4 m and x c = 5 m. Figure 1 
plots the pressure coefficient distribution for this c a se . 
The strength of the pressure gradient is measured by the 
nondimensional pressure gradient parameter 


ft 5 * dp 


(26) 


where 5* is the displacement thickness. An average 
value of (Ij. = 20 was maintained in the equilibrium 
region. 

A zero pressure gradient flat plate boundary layer is 
used to match Re ^ , the momentum thickness Reynolds 
number, in the test section. Figure 2 plots the computed 
and experimental Reynolds numbers based on the length 
scales. Good agreement is maintained with the integral 
length scales. 

In order to mimic the conditions stipulated by Strat- 
ford, the added diffusion term was implemented in the 
boundary layer code as 

X«(X) (27) 


where 



H (x) is the Heaviside function, i.e. H(yJ = 1 for % > 0 
and zero otherwise. 

Skin friction values are compared in Figure 3. Skin 
friction is defined as 
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C ft = (29) 

-p U 2 
2^ ee 

where ( ) c denotes a quantity based on boundary layer 
edge values. k -(0 results are about 50% larger than the 
expected results. The new model matches the data quite 
well with only a slight underprediction in the equilib- 
rium pressure region of approximately 10%. 

The mean velocity profiles are compared in Figures 
4 and 5 with the law of the wall, where 



(30) 


where v denotes the kinematic viscosity and the friction 
velocity is defined as 



The data and the new model results match quite well, 
whereas the k ~ CD model underpredicts the profiles. The 
log law maintains its validity in adverse pressure gradi- 
ent flows for a range of y + . The pressure gradient affects 
the point of departure from the log law. 

Figure 6 plots the shear stress in the boundary layer 
at all six measurement stations. The peak of the shear 
stress is at approximately y/8 = 0.5 . In a zero pressure 
gradient flow, the maximum shear stress is at the wall. 
Similar behavior was found in other adverse pressure 
gradient flow casesA®* 12 The modified k - co results 
match the data quite well for y/5 ^ 0.5 and only 
slightly overpredicts the outer region. Both methods find 
the correct location of the peak, but the magnitude is 
consistently underpredicted by the k < o model. 

Based on these results, we notice that z w , as deter- 
mined by the original fc-co model, is gready over-esti- 
mated in regions of strong adverse pressure gradients. In 
such regions, the velocity distribution described by 
Stratford for vanishing shear stress adequately describes 
the near wall behavior and thus matches the skin fric- 
tion. 


Experiment (1\ 

Experiment (2) consists of data from an adverse 
pressure gradient, supersonic turbulent flow on a flat 
plate. At the start of the test section, the reference Mach 
number is M re f = 2.92. A fit of the wall pressure is 
shown in Figure 7. The maximum j3j = 5.8 in this pres- 
sure field, much lower than the values found in Experi- 
ment (I). As before, a zero pressure gradient flat plate 


boundary layer is used to match the initial conditions. 
Emphasis was placed on matching Re at the start of 
the of the pressure gradient. The error for Re 5 „ , Re h , 
and Re B are 0%, 2%, and 17% at the start of the test 
section. 

The skin friction predictions from both models are 
compared to the data in Figure 8. The skin friction is 
defined here as 


'/.ref 


2 P relief 


(32) 


where ( ) r< ^ denotes the reference values at x c - L The 
wall shear stress is determined by the authors based on 
the Bradshaw and Unsworth 13 calibration of a Preston 
probe, and also by solving for the friction velocity based 
on either a Van Driest 14 or a Carvin ex ai 15 velocity 
transformation. The three methods agree within 6% of 
each other. The £-0) model outperforms the new model 
and matches the data well. The new model underpre- 
dicts the skin friction by almost 22%. 

The typical mean velocity profile is plotted in Fig- 
ure 9. The k ’ co model matches the profile in the near 
wall region better than the present model. However, in 
the outer region, both methods under-predict the slope 
for the velocity rise. The mean density profile is com- 
pared with the theory in Figure 10 for the same location. 
Both models overestimate the density increase in the 
near wall region and the final density value at the edge 
of the boundary layer. 

Figure 1 1 displays the calculated and experimental 
shear stress for two locations. A dramatic difference 
between the two models is seen at x c = 1 .254. The stan- 
dard k-t o model follows the data trend in the outer 
region well, but overshoots the data peak by approxi- 
mately 18%, which is within -5% to 30% uncertainty 
associated with these data. Also note that for y/8 < 0.5 , 
the data is considered problematic. The current model 
misses the proper trend from the wall out to y/8 = 0.7 
with a -12% difference in the peak value. 

The final profile point at x c — 1.381 shows better 
agreement between the new model and the data in the 
outer region, but once again the trend of the dafa is not 
well followed. The k - co model performs better than the 
other model in locating the position of the maximum 
shear stress. The error bounds mentioned for the previ- 
ous profile point apply here as well. Thus, the k - co model 
prediction falls relatively close to the experiment. 

These results find the £-co model superior to the 
new formulation at supersonic speeds, as well as to the 
fc-e model applied to this same data by Becht and 
Knight 16 . 

The inability of the present model to predict super- 
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sonic boundary layer behavior is expected, since 
adverse pressure gradients affect subsonic and super- 
sonic boundary layers differently. Although the k- a) 
model performed well for Experiment (2), it is well 
known that none of the existing models perform well in 
regions of shock induced separation. As a result, a dif- 
ferent diffusion model is necessary to improve existing 
predictions. This will require additional well docu- 
mented experiments. 

Concluding Remapfo 

The major contribution of the present work lies in 
the fact that it identifies inadequate diffusion modeling 
as the cause for the inability to model separated flows or 
flows on the verge of separation. The results suggest that 
a new diffusion mechanism needs to be introduced in 
the k - equation in order to explain observed behavior of 
turbulent flows in the presence of strong adverse pres- 
sure gradients at low speeds. The current modification 
of the £-(o model with an additional diffusion term based 
on the pressure gradient showed good agreement with a 
strong adverse pressure gradient in an incompressible 
flow. However, because subsonic and supersonic bound- 
ary layers behave differently in the presence of adverse 
pressure gradients, a more complete diffusion model is 
needed. This will require additional measurements and 
additional modeling. 
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Figures 



Figure 1 . Pressure distribution and curve fit for Expert- 
ment (1). 



Figure 2. Comparison of length scale Reynolds num- 
bers for Experiment ( 1 ). 



Figure 3. Comparison of skin friction for 
Experiment (1), 



Figure 4. Comparison between £-co model and data of 
the mean velocity profiles for Experiment (1). 



Figure 5. Comparison between present model and H^ ta 
of the mean velocity profiles for Experiment (1) 
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Figure 6, Comparison of shear stress profiles for Experiment (1). 
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Figure 7. Pressure distribution and curve fit for Experi- Figure 10. Comparison of mean density profile at x c 
ment (2). 1 ,254 m for Experiment (2). 



Figure 8. Comparison of skin friction values for Experi- 
ment (2). 





Figure 11. Comparison of shear stress profiles for 
Experiment (2). 


Figure 9. Comparison of mean velocity profile at x c = 
1.254 m for Experiment (2). 
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A k - C (ENSTROPHY) COMPRESSIBLE TURBULENCE MODEL 
FOR MIXING LAYERS AND WALL BOUNDED FLOWS 

G. A. Alexopoulos *, and H. A. Hassan f 
North Carolina State University, Raleigh, North Carolina 27695-7910 


Abstract 


A compressible turbulence model based on the k - £ 
(enstrophy) model of Robinson et al. is developed to cal- 
culate free shear and wall bounded flows. The model is 
validated by comparing to the mixing layer data compiled 
by Settles and Dodson. It is shown that it reproduces 
measured similarity profiles of velocity, turbulent kinetic 
energy and shear stress. Moreover, its predictions are su- 
perior to the k - € model with traditional compressibility 
corrections. Furthermore, the model predicts skin fric- 
tion, heat transfer rates and law of the wall profiles with 
good accuracy for a variety of wall bounded flow cases. 


Introduction 

Currently, the two most used two-equation turbulence 
models are the k - e and the k - u models. Both of these 
models are not successful in predicting free shear flows 
and wall bounded flows with the same set of model con- 
stants and boundary conditions* . Part of the difficulty is 
a result of the fact that such models were originally de- 
rived for flows characterized by high turbulent Reynolds 
numbers. Although no adjustments are made in the tradi- 
tional k - u> model at low turbulent Reynolds numbers* , 
the adjustments made in the k - e models are usually 
guided by asymptotic expansions and not by the exact 
form of the e-equation. Because of this, it is believed that 
available low turbulent Reynolds number corrections em- 
ployed in the k - e model are not consistent with the exact 
e-equation and thus lack the correct physics. 

The above considerations were the basis of a new two- 
equation k - £ (enstrophy) turbulence model developed 
by Robinson et al. 2 for incompressible flows. Instead 
of modeling the exact equation for dissipation, attention 
was focused on modeling the exact equation for the vari- 
ance of vorticity, or enstrophy. Thus terms that depend 
on the mean flow, which are negligible at high turbulent 
Reynolds numbers but not a low turbulent Reynolds num- 
bers, were retained and modeled. It is shown in Ref. [2] 
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that the k - ( model reproduces all available measure- 
ments of growth rates and turbulent shear stress distri- 
butions for a variety of free shear layers using the same 
set of model constants. 

The objective of this work is to extended the k - C model 
of Robinson et al. 2 to compressible turbulent flows. The 
success of the k - £ for reproducing a variety of incom- 
pressible flows should provide a good basis for the cur- 
rent model. As in the incompressible model, the current 
model will be based on the exact compressible equations 
to ensure that the correct physics is incorporated. 

Compressible flows require for their description a ve- 
locity field and two thermodynamic variables such as the 
density and temperature. Because of this, the fluctua- 
tions of these thermodynamic variables are as important 
as those of the velocity in determining the resulting tur- 
bulent flow. As a result, traditional two-equation and 
stress models 3 have proven to be inadequate in describ- 
ing such flows. Therefore, it appears that an appropriate 
compressible turbulent flow model of the “two-equation” 
variety, should include six equations that describe vari- 
ances of velocity, density and temperature and their re- 
spective dissipation rates. 

A simplification of the above approach can be achieved 
by invoking Morkovin’s hypothesis 4 . According to this 
hypothesis, the pressure and total temperature fluctua- 
tions are small for non-hypersonic compressible boundary 
layers with conventional rates of heat transfer, i.e. 


As a result, 





(7 — 


(1) 

( 2 ) 


In the above equation, p } T' 0 ,p , 7 1 ' and u are the fluc- 
tuating pressure, total temperature, density, temperature 
and velocity. Also M is the Mach number and 7 is the 
ratio of specific heats. The remaining variables represent 
mean properties. Based on Eqn. (2), equations governing 
variances of turbulent quantities can be taken as the equa- 
tion for the turbulent kinetic energy. However, equations 
governing the dissipation rates of the resulting variances 
may not be the same. 

The approach employed in developing a k - £ compress- 
ible model is to use the guidelines set in Ref. [2] to model 


1 



the additional terms. Thus, for low Mach numbers the 
new model reduces to that of Ref. [2]. It should be noted 
that existing k - t models assume the effects of compress- 
ibility to be limited to the ^-equation. This is not the 
case because the exact compressible e-equation is differ- 
ent from that for incompressible flows. 


Governing Equations 

The mean-flow equations governing the conservation of 
mass, momentum and energy can be found in Ref. [1]. 
The compressible turbulent kinetic energy equation 1 
(Fav re- aver aged) and enstrophy equation can be written 
as 
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with being the eddy (turbulent) viscosity. 


where 
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c = w-w'. 
p 
p 




J \) 
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Here c is the turbulence energy dissipation rate, p is the 
viscosity, Q, and u i are the mean and fluctuating vortic- 
ity, tij and r t; - are the laminar and turbulent (Reynolds) 
stress and e IJ/t is the permutation tensor. 

The incompressible terms of Eqns. (3) and (4) were 
modeled in Ref. [2] with the remaining compressible terms 
modeled in the Appendices A and B. Upon modeling, 
Eqns. (3) and (4) take the form 


Dk_ _ dUi_ d 
P Dt - T ,j dxj + dx. j 


(d th\ dk_ 
\ 3 + a k ) d Xj 


Results and Discussion 

To validate the present model, comparisons were made 
with various experimental measurements of compress- 
ible mixing layers and wall bounded flows. The model 
constants for the incompressible terms are tabulated in 
Ref. [2] and these were unchanged. The new terms are 
modeled in Appendices A and B. For the flows consid- 
ered here, mixing layers and flat plates, only two terms 
were found to be significant. The first term appears in 
the ^-equation (see Appendix A) and represents dilata- 
tional dissipation. The second term is in the (^-equation 
(see Appendix B), which is a production term 

resulting from dilatational effects. The remaining terms, 
which may be important for other complex flows, are in- 
cluded for completeness. The model constants for the new 
compressible terms which were determined by numerical 
optimization are 


Ci = 0.6 , C Cl = 1.35 

The ability of the model to calculate compressible mix- 
ing layers will be illustrated first. The two mixing layer 
cases that were chosen survived the scrutiny of Settles and 
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Dodson 0 . The first set is due to Samimy and Elliott 0 and 
Elliott and Samimy 7 while the second set is that of Goebel 
and Dutton 8 . Settles and Dodson tabulate and plot (see 
Table II and Fig. 1 of Ref. [5]) normalized mixing layer 
growth rates versus convective Mach numbers. The data 
is collected from various sources. As may be seen from 
the data and figure, there is a great deal of scatter render- 
ing such tabulations and plots meaningless. The reason 
for the scatter can be traced to a lack of uniformity in 
calculating and defining growth rates. For most of the 
available data, the growth rate db/dx ^ b/x, where b is 
the width of the mixing layer. Moreover, there is a lack 
of uniformity in defining the width of the mixing layer. 
Because of this, judging the worth of a compressible tur- 
bulence theory based on its prediction of growth rates is 
somewhat misleading. 

Rather than restrict our comparisons to growth rates, 
we opted to compare the predictions of the theory with 
the measured velocity, turbulent kinetic energy and tur- 
bulent shear stress. Since the k - e model performs bet- 
ter than the k - w model for free shear flows 1 , compar- 
isons were made with a k - e model using the compress- 
ibility correction model proposed by Wilcox 1 . It may 
be recalled that this model incorporates the best feature 
of both Sarkar et al. 9 and Zeman's 10 models. Table 1 
summarizes the cases given in Ref. [5] and described in 
Refs. [6] - [8]. 

Comparison of present theory with the k - e model and 
experiment are shown in Figures 1-7. The computational 
procedure is similar to that employed in Ref. [2]. Taking 
compressibility into consideration, the similarity form of 
the variables for mixing layers is 

f pyty 

^ ~ P\x 

v{x,y)z=e&. [V(»7) - | a*W(i»)] 

*(*,y) - U?K(r,) 


u(x,y) = UiU(ri) 
h(x,y) = h\H(y) 
C (x,y) = ^E(rj) 


Y — r* KlLn) 
Iy “ E(n) 


where the symbols have the same meaning as in Ref. [2]. 

In presenting the data, we followed the suggestion of 
Ref. [6] and scaled y with <5^ , the vorticity thickness, de- 
fined as 


bw — 


A U 



A U ~ U x - U 2 


( 9 ) 


where subscripts 1 and 2 represent the fast and slow 
streams, respectively. Each figure consists of three plots 
for normalized velocity, turbulent kinetic energy and tur- 
bulent shear stress, r xy . 

It is seen from the figures that the present theory is 
in good agreement all cases of Refs. [6] and [7] and with 
cases 2, 3 and 4 of Ref [8]. For these cases, the k - e model 
shows poor agreement with experiment. The agreement 


is not good with case 5 of Ref. [8], which is at the highest 
stagnation temperature ratio and convective Mach num- 
ber. The reasons for this are not clear at this time. Note 
that this case is in good agreement with the k - t model. 
Furthermore, applying the compressibility correction of 
Wilcox to the k - £ model yields worse agreement with 
experiment than the k - e with compressibility corrections. 

To illustrate the ability of the model to compute more 
complex flows, three wall bounded flow cases were cho- 
sen and are summarized in Table 2. The first two cases 
are adiabatic wall cases obtained from the compressible 
turbulent boundary layer data compiled by Fernholz and 
Finley 11 . The first case is from from Shutts et al. 12 at 
Moo = 2.24 in air while the second set of data is that of 
Watson et al. 13 at = 10.31 in helium. The third set 
of data is a cold wall case at = 8.18 in air extracted 
from Kussoy and Horstman 14 . These cases were a subject 
of a previous investigation by Sommer et al. 15 

Calculation of the wall bounded flows was obtained us- 
ing a marching boundary layer code, where the model 
boundary conditions are k^ = k w = ^ = 0 and 
|ty= 0. The eddy viscosity is defined as 

pit 2 

= , Cp = 0.09 (10) 

where 

U = min(/;, 1.0) 

r '’{ l + < n > 

and 

y + = ^ , u T = , C lii = 0.5 , = 29.5 (12) 

Comparisons of the k - £ model with the experiment 
are shown in Figures 8 - 10. In presenting the results, we 
followed the suggestions in the critical review of the data 
by Fernholz and Finley 16 . It was determined in the re- 
view that certain differences were observed in the data 
that were due to a sensitivity to the plotting method 
rather that from an observed physical phenomenon. Fur- 
thermore, it was suggested that a plot of T/T& versus 
U/Ub represents a valid comparison for temperature for 
examining the Crocco - van Driest relation. The first two 
experimental cases use the Crocco - van Driest relation 
by assuming constant total enthalpy and pressure across 
the boundary layer to determine the temperature profiles 
from the measured velocity profiles. As a consequence, 
the temperature data is not an independent measurement 
for these two cases. 

Each figure consists of a plot of temperature versus 
velocity distribution and a law of the wall comparison 
between the experiment and present theory. The calcu- 
lations for each case were performed by matching the re- 
spective experimental fte*, given in Table 2. It can be 
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seen from Figures 8 and 9 that the k - Q model predicts 
the experiment with good accuracy. However, there is 
a slight discrepancy between theory and experiment for 
the law of the wall distribution of Figure 9. It was de- 
termined in Ref. [16] that this case may not have been 
fully turbulent, which may explain the differences in the 
law of the wall. The difference in measured Cj (given in 
Table 2) and theory for these two adiabatic wall cases is 
0.1% and 1.5%, respectively. 

In Ref. [14] the measurements of velocity and temper- 
ature are obtained independently. Figure 10 represents 
the cold wall case comparisons and is showing an excel- 
lent agreement between experiment and present theory. 
Further, the difference in measured Cj and Ch for this 
case is 5.6% and 8.6%, respectively, which is well within 
experimental error. 

Conclusions 

Comparison between experiment and theory show good 
agreement in the calculation of compressible mixing lay- 
ers. Furthermore, the current k - £ was shown to be 
generally superior to k - e model in predicting mixing 
layer experiments. Finally, an examination of compress- 
ible wall bounded flows show an excellent agreement in 
the calculation of skin friction and heat transfer rates and 
a good prediction of temperature profiles and law of the 
wall distribution. Future work will involve further valida- 
tion by implementing the model in a Navier-Stokes solver 
and calculating more complex flows. 
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Appendix B: Modeling of the 
C Equation 


Appendix A: Modeling of the 
k Equation 

The dissipation rate in Eqn. (3) is defined as 


1 ' ' 

1 . ui x S n 

The above first order tensor quantity has units of [T] 
and is modeled as 


pc = v 


+ 2pu- j u- i - -pu'-u'l- 


(A.l) 


> du 'j ~ M~ 

= n, 


(B.l) 


"2 . ^ / // Si 

P^i +3 P(«i,i) 


where l/r p was defined in Appendix A and Q = 


+ 2 i/ 
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The second term in Eq. (A.2) can be written as 
d 


The above scalar quantity has units of [jj] and is mod- 
eled as 


2i/ 


dxi 


d / „ ~ \ ..du'j 

jj - (?u, ) - 2 P u, 


'2 ' .M ^ U j __ ^ 

‘ Ji ‘ dx 3 ~ C ^2 t„ 
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~ 2 


0 

0 *i 


f d ( — TT\ 

„ du" } 


""•4/ 


(A.3) 


and is normally neglected, consistent with the assump- 
tions of homogeneous turbulence. In this work, the term 
is included in the diffusion term. 

The approach employed in modeling the various terms 
follows closely that employed in Ref. [2]. Moreover, model 
constants developed in Ref. [2] were unchanged in this 
model. 

The compressibility term in the k equation has gener- 
ated a great deal of discussion lately and led to the devel- 
opments of models by Sarkar et al, 9 , Zeman 10 , Wilcox 1 
and Ristorcelli 17 among others. None of these models 
were adopted here because they do not blend smoothly 
with the incompressible limit. Moreover, there is no rea- 
son to assume that the time scale governing this term is 
dependent on Instead the term 


q 'K 

The quantity t km is defined as 

/ / du k dx i* \ du : 

^ = + < B - 3 ) 


where A = — |/i. Thus, 


d 2 U; 


km _ ( & u l fc . d U m \ m a 

p \ 'oZT' + o ' :: ■ _ I + A b krn 7r— 1 ~ ( B - 4 ) 

dl m Xj 


dx m y dx'l dif-x 
'd 2 u‘ . dV 


= p 
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d 2 U; 


dx 2 m ' dx k x m ) +X dx k ij 


v ^p( u i.i) 2 

(A.4) 

k 

Ci — 

Tp 
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and 


is modeled as 


where C\ is a model constant and r p is a time scale. Be- 
cause the term under consideration is zero when the den- 
sity is constant, t p is modeled as 


. dt' k , d 2 u‘ k 

UjkWi d^Z = (l€iikUi d^ 


, d 2 u '■ 

+ (\ + fi) Cijtuii Q x (B-7) 


= 2 /ir; 
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where 2 r k - _ eij k w i . Using the gradient diffusion approx- 
imation, the above terms can be modeled as 


o _ n fh dc 

P kl dxl 4 


2 dx 3 

,,, , v . C Vt dp 

Hence, the term in question is modeled as 


fuJ dt ' km -n, ^ 5( > O- to dp 

,;t <'^7 -Q4 “~ + a5 ~ 
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2 dx. 
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Since 




,dp_ 

dxj 


is a vector, it is modeled as 


= —2 r 




Ujk^i 


, dp 
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os dp 
Tp dx k 


(B. 12) 


(B. 13) 


Table 1: Mixing Layer Experiments 


Case 

Ui/Ui 

pi!p\ 

^OJ /To 2 

To,,r. a (K) 

M \ , M 2 

M c 

Samimy and Elliott - Case 1 

0.36 

0.64 

1.00 

276, 276 

1.80, 0.51 

0.51 

Samimy and Elliott - Case 2 

0.25 

0.58 

1.00 

276, 276 

1.96, 0.37 

0.64 

Elliott and Samimy - Case 3 

0.25 

0.37 

1.00 

276, 276 

3.03, 0.45 

0.86 

Goebel and Dutton - Case 2 

0.57 

1.55 

1.96 

578, 295 

1.91, 1.36 

0.46 

Goebel and Dutton - Case 3 

0.18 

0.57 

1.00 

285, 285 

1.96, 0.27 

0.69 

Goebel and Dutton - Case 4 

0.16 

0.60 

1.24 

360, 290 

2.35, 0.30 

0.86 

Goebel and Dutton - Case 5 

0.16 

1.14 

2.25 

675, 300 

2.27, 0.38 

0.99 


M c — (f/i U c )/ a i , U c — { a \U2 o.2Ui)/(d\ 0 . 2 ) 


Table 2: Wall Bounded Experiments 


Case 



(K) 

Res 

T w /T aw 

C/xlO -3 

C\xl0 3 

Shutts et al. 

2.24 

0.2364 

339 

V 


1.62 


Watson et al. 

10.31 

7.5576 


■MIMm 

1.0 

0.24 



8.18 

6.0795 

1166 


0.3 

0.98 
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Figure 3: Turbulence Model Comparisons with Experi- Figure 4 : Turbulence Model Comparisons with Experi- 
ment (Elliott and Samimy Case 3) ment (Goebel and Dutton Case 2) 
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Figure 7: Turbulence Model Comparisons with Experi- 
ment (Goebel and Dutton Case 5) 
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Figure 8: Comparison of Calculated Mean Temperature 
- Velocity Distribution and Law of the Wall versus Exper- 
iment (Shutts et al.) 





Figure 9: Comparison of Calculated Mean Temperature Figure 10: Comparison of Calculated Mean Tempera- 
- Velocity Distribution and Law of the Wall versus Exper- ture - Velocity Distribution and Law of the Wall versus 
iment (Watson et al.) Experiment (Kussoy and Horstman) 
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Abstract 

The two-equation turbulence model of Robinson et al. 
is examined for a variety of wall bounded cases including 
a boundary layer solution of a flat plate and a Navier- 
Stokes solution of airfoils. It is shown that traditional 
methods of modeling turbulent diffusion are not appro- 
priate for describing turbulent separated flows. A ten- 
tative proposal for modeling such flows is presented and 
the predictions of the model are compared with measure- 
ments for a difficult case involving an NACA 0012 airfoil. 
The results presented here clearly illustrate that a two- 
equation turbulence model can be developed to describe 
a variety of flows, for which a two-equation model is ap- 
propriate, using the same set of model constants. 

Introduction 

Existing k-w and k-e turbulence closure models are in- 
capable of describing wall bounded and free shear flows 
using the same set of model constants and boundary con- 
ditions. It is generally agreed that the problem with ex- 
isting models comes from the dissipation equation, c or w. 
The reason for this can be traced to the fact that these 
equations were developed for high turbulent Reynolds 
numbers, Re t} but are being employed in situations where 
Re t is low. 

For most turbulent shear flows Ret is typically large 
enough so that one can assume the small scales are nearly 
independent of the large scales and thus their dissipation 
rate is isotropic. If one makes this isotropic assumption 
then the terms in the exact dissipation equation which 
depend on the mean flow can be neglected. However, for 
most shear flows, the turbulent Reynolds number varies 
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from very large values in the bulk of the flow to very small 
values near walls and outside of the boundary layer. 

For regions of low Re t , the small scales of turbulence 
are weakly dependent on the large scales, and therefore, 
on the mean flow. Because of this, terms appearing in 
the exact dissipation equation which depend on the mean 
flow cannot be neglected. 

The above considerations were the basis of a new two- 
equation model 1 . Instead of modeling the exact equation 
for dissipation, attention was focused on the enstrophy or 
the variance of vorticity equation 2 . The model developed 
in Ref. [1] was used to describe free shear layers (wakes, 
jets, and mixing layers). The same set of model con- 
stants was used in all calculations. Excellent predictions 
of growth rates, shear stress, and velocity distributions 
were obtained. The object of this work is to use the same 
set of model constants developed in Ref. [1] to study wall 
bounded flows. The implementation of the model will 
proceed in a number of stages. In the first, the model 
is implemented in a boundary layer code 3 and the re- 
sults are illustrated by calculating a flow past a flat plate 
and its wake. Second, the model is implemented in a 
Navier-Stokes code 4 and the results are illustrated by 
calculating attached flows past airfoils. In the third, at- 
tention was devoted to separated flows past airfoils. It 
was determined that traditional models of turbulent dif- 
fusion are not adequate for calculating separated flows. 
This resulted in the development of a tentative diffusion 
model. 

The model was validated by calculating flows past an 
NACA 0012 at an angle of attack, a, of 2.26 deg and 
a free stream Mach number, Afoo, of 0.79. It may be 
recalled 5 that existing two-equation models were unable 
to predict the pressure distribution for this case; only the 
Johnson and King model 6 performed well for this case. 
The results presented are intended to demonstrate that 
one set of model constants can be used to describe all 
types of shear flows. Obviously, applications were limited 
to those cases where a two-equation model is expected to 
give good results. 
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Analysis 

Assuming incompressible flow, the governing equations 
are the Reynolds averaged continuity, momentum, turbu- 
lent kinetic energy, and enstrophy. These equations can 
be written as 

dUi n 

= ° <» 
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6 vyK P 

with p t being the eddy viscosity. The model constants 
are provided in Table 1. 

Two departures from earlier modeling should be noted. 
First, when testing the model in a Navier-Stokes code, it 
was determined that the first term in equation 7 should 
be replaced by 

[ divtCij) d(v t n,y 

(Jr dxj dxj dx { 


(9) 

The original modeling was consistent with Taylor’s vor- 
ticity transport theory However, we encountered simi- 
lar difficulties here. The above modification removes the 
difficulty. It should be noted that, because of similarity 
considerations, the term under discussion did not con- 
tribute to the similarity equations employed in Ref. [1]. 
Therefore, the above modification does not change any of 
the model constants. 

The second adjustment concerns the dissipation term 
in the ( equation. The turbulent time scale ~ can not 
be less than Kolmogorov’s time scale. To allow for this 
result, the dissipation term now has the form 


Rk + $ 

For the results presented here 6 = 0.1. 


( 10 ) 


p is the density, p is the viscosity, u-, u;-, and p' are 
the fluctuating velocity, vorticity, and pressure, and 17, 
is the mean velocity. Equations 3 and 4 were modeled 
in Ref. [1] and the results can be written as 



Results and Discussion 
I. Attached Flows 

For high turbulent Reynolds numbers, the eddy viscos- 
ity is chosen as 


C k 2 

= , C/j = 0.09 (11) 

For wall bounded flows the eddy viscosity and the tur- 
bulent Reynolds number are both zero at the wall. Be- 
cause of this, the above expression is traditionally multi- 
plied by the function /^. There is no unanimity on the 
form of in the near wall region 8 . More recent im- 
plementations such as that of Speziale et al 9 employ a 
function that depends on both Ret and where 
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where r w is the wall shearing stress. Because separation 
may take place in the presence of adverse pressure gradi- 
ents a different representation that does not require y* is 
employed. The resulting expression is given as 

/„ = min(/;,1.0), /; = (n. £k) exp (.3, 


the expression for the velocity in the log- law region, and 
the manner in which k varies with y in the near wall re- 
gion [See Table 4.5 in Ref. [10]]. Recall that the correct 
near wall variation of k can be enforced through the use 
of the second boundary condition on £. 

Figures 1- 3 show plots of Jb+, £+, and (-puV) + vs. 
y + in the near wall region where 


where 


C Ml =4.0 , 40.0 (14) 

The boundary conditions for the two-equation model 
are 


k w = 0 

koo = ^(TW ro ) 2 (15) 

and 

<9C , 

U = 0 , or (16) 

d (I8k\ 

= c,u aty=0 

where T is the turbulent intensity (« 1%) and Coo is de- 
termined by specifying a free stream ratio of ^ along 
with too. The first boundary condition for Cw amounts 
tto a simple extrapolation when a first order difference for 
In, is used. The second boundary condition for C is ob- 
tained by examining equation 6 in the near wall region. 
This equation is then used to solve for Cu# as follows. First, 
we expand as a central difference 


c«,= 


dk 

_ fly 


i dk 
I 1 dy 

3A y 


where T’ represents the first point off of the wall and ’w* 
represents the wall. Since k ~ y 2 near the wail, we have 


dk 

dy 


U=o 


(17) 


and, therefore, 


3 (Ay) 2 3 (Ay) 2 


(18) 


Using this boundary condition for £ forces a )b ~ y^ varia- 
tion near the wall. Results for both boundary conditions 
are presented for comparison. 

There are a number of tests a turbulence model must 
meet to be considered a successful model. First, and fore- 
most, it must predict the correct skin friction and pres- 
sure coefficients and other near wall measurements ® of 
shearing stress, kinetic energy, and dissipation. Second, it 
must predict the correct growth rate, asymptotic shearing 
stress distribution, and velocity in the far wake. Finally, 
it must predict the correct J B’ constant that appears in 


k+ 



—pu'v* 


Ihl 

U? 


(19) 


Also shown is the ” average” of available data 8 together 
with the experimental measurements. The predictions 
compare well with the results summarized in Ref. [8] and 
experiments. It appears from Figure 1 that the peak *+ 
is under predicted, nevertheless, it lies within the range 
of available experimental data (Jt+ « 2.8 - 6.0). Figure 4 
compares the calculated skin friction coefficient with the 
correlations of Cole and Schoenherr 11 . Again good agree- 
ment is indicated. 

The constant ’B* in the log-law correlation 


«+ = i/ny+ + B u i (20) 

is typically 5. A plot of(uC-l In y+ shows that , B’ is not 
exactly constant, but rather, varies slightly from approx- 
imately 4.4 to 5.2 over a y + ranging from 50 to 500. 

The next set of figures compares calculated and mea- 
sured wake data. Figures 5- 7 show a comparison of nu- 
merical, asymptotic, and experimental results for the de- 
fect velocity at the centerline, & [U e is the edge veloc- 
ity, W 0 is the centerline defect velocity], wake half width, 

b, and the maximum shear stress (ppr), where r is the 
shear stress, as a function of f where 6 is the momentum 
thickness. The asymptotic solution is that of the far wake 
with a constant eddy viscosity. The experimental results 
are obtained from Pot 12 and Weygandt and Mehta 13 . 
Again, good agreement is indicated. 

Figures 8- 9 compare the defect velocity, W, and the 
shear stress distributions across the wake. Here, also good 
agreement with experiment 12 ’ 14 is indicated. 

The remaining figures show the Navier-Stokes (NS) so- 
lution for a NACA 0012 airfoil. All of the NS solutions 
were run using a third order upwind biased Roe scheme. 
Figure 10 shows the grid (301x101) which has an initial 
normal spacing of l.e~ 6 C and an outer boundary of 15C\ 
Figures 11- 14 show the pressure and skin friction coeffi- 
cients compared with the experimental data 15 and with 
the standard k-u model of Wilcox 10 . The Wilcox model 
was chosen because it is widely known that this model 
gives good agreement with experiment for wail bounded 
flows. 

Figures 11- 12 correspond to a Mach number of 0.502 
and an angle of attack of 2.06 deg. Once again, good 
agreement for both C p and Cj e are shown. For this case, 
the choice of boundary condition on £ affects the nu- 
merical stability of the solution by reducing oscillations 
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around the trailing edge. Except for this small region 
around the trailing edge, the choice of boundary condition 
has very little effect on the solution. This insensitivity 
was also noted when examining the flat plate boundary 
layer solution. The reason for this insensitivity is seen by 
noting that for a flat plate, the k-£ model (with simple 
extrapolation for £), produces a k ~ y 2 behavior near the 
wall. Therefore, using the second boundary condition to 
impose a y 2 variation should have no effect. 

The next case (Mach=0.80, a = 0.0) examines the be- 
havior of the k-£ turbulence model for a transonic airfoil. 
From figure 14, it is noted that an extrapolation bound- 
ary condition for £ causes a large spike in the skin friction. 
This oscillation is eliminated by imposing the second be 
for Therefore, it appears that bc2 is the appropriate 
choice for a boundary condition for non-simple flow fields. 
When using bc2, good agreement with the k-w model is 
shown. This boundary condition was used in calculating 
separated flows. 


II. Separated Flows 

It is a known fact that both k-w and k-t models are, 
in general, unable to predict separated flows 16 . In a 
recent investigation, Rao and Hassan 16 showed that the 
reason for this behavior is a result of inadequate modeling 
of the diffusion term in the k-equation. There are two 
ways of modifying the diffusion term: one by remodeling 
it, which is the approach followed in Ref. [16]; the other 
by adjusting the model ’constant’ in the diffusion term. 
Because subsonic and supersonic boundary layers behave 
differently in the presence of adverse pressure gradients 
it has not been possible to remodel the diffusion term 
in such a way so as to accommodate both subsonic and 
supersonic flows. As a result, the second approach, which 
entails adjusting the model constant, was employed here. 

The approach presented here is tentative because it 
does not have the desired general form. It entails adjust- 
^ in the turbulent kinetic energy equation by sub- 
tracting from it a dimensionless pressure gradient param- 
eter. A number of guidelines were set for selecting such a 
parameter. It may be recalled that boundary layer theory 
is rich with such parameters, unfortunately, most of them 
are not convenient to compute in a Navier-Stokes solver. 
Moreover, it is desirable to avoid quantities that require 
a search in the flow field. Further, the new model should 
be such that it has little or no effect on the attached flow 
cases discussed above. Finally, the model must be able to 
handle a difficult separated flow situation like the NACA 
0012 at a = 2.26 and MV* = 0.8. 

In the proposed model, in the k-equation is replaced 
by 

“A ( 21 ) 


where 

^ = qZdS ’ q °° ~ 2 PooU °° ( 22 ) 

and 

0 < 0t < 1*4 (23) 

It was necessary to limit (3 t because the term has a 
large spike around the shock. Note that when = 0, 
the original model constants are recovered. 

The above model was used to calculate the flow past an 
NACA 0012 at a = 2.26 and Mqq = 0.79. The results are 
shown in Figures 15 and 16, and compared with a calcu- 
lation by Rumsey and Anderson 17 using the Johnson and 
King model, and a k-w model. Figure 15 compares the 
pressure distribution while Figure 16 compares the skin 
friction. The present and Johnson and King give good 
agreement with experiment, but the k-w does not. Un- 
fortunately, a lack of experimental data does not permit 
a meaningful discussion of skin friction results. 

Conclusions 

From the results indicated here, it has been shown that 
the model developed in Ref. [1] can be used for both wall 
bounded and free shear flows using one set of model con- 
stants. It has also been shown that the model can be eas- 
ily incorporated into either a boundary layer or Navier- 
Stokes code. Finally, it has been demonstrated that the 
second boundary condition for ens trophy is more accu- 
rate for non-simple flow fields. We recognize that, al- 
though the tentative diffusion model presented here is a 
good start, a more general model is required. Moreover, 
additional validations would be desirable especially for 
separated flows. 
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Table 1: Turbulence Model Closure Coefficients 


Constants 

*-< 


0.09 

K 

0.40 

<*3 

0.35 

1 P* 

0.42 

06 

2.37 

06 

0.10 

07 

ro.75 

06 

2.30 

<Tr 

2.00 

i 

<7fc 

1.80 


1.46 
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Figure 1: Near Wall Behavior of Turbulent Kinetic Energy (Flat Plate) 



Figure 2: Near Wall Behavior of Enstrophy (Flat Plate) 



Figure 3: Near Wall Behavior of Reynolds Stress (Flat Plate) 
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k-t (E/isroony) 

□ Col* (1956) 

O Scho*nh*rr (1932) 
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Figure 13: Pressure Distribution for NACA 0012 Airfoil (M» = 0.8, Re = 9.0e6, a = 0.0) 



Figure 14: Skin Friction Comparison for NACA 0012 Airfoil ( M ^ = 0.8, Re ^ = 9.0e6, a = 0.0) 



